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Hyperfine anomaly in Be isotopes in the cluster model and the neutron spatial 

distribution. 
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The study of the hyperfine anomaly of neutron rich nuclei, in particular, neutron halo nuclei, 
can give a very specific and unique way to measure their neutron distribution and confirm a halo 
structure. The hyperfine structure anomaly in Be + ions is calculated with a realistic electronic 
wave function, obtained as a solution of the Dirac equation. In the calculations, the Coulomb 
potential modified by the charge distribution of the clustered nucleus and three electrons in the 
ls 2 2s configuration is used. The nuclear wave function is obtained in the core+nucleon model of 
9,11 Be. The aim of this study is to test whether the hyperfine structure anomaly reflects a halo 
structure in 11 Be. 
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I. INTRODUCTION 

Exotic (halo) nuclei are currently a subject of inten- 
sive experimental and theoretical studies. The interest 
in neutron-rich nuclei is to a large extent driven by ex- 
perimental facilities and new experimental methods for 
studying the nuclear matter distribution. In particular, 
the ion trap method pi and the NMR (Nuclear Magnetic 
Resonance) methods j2j allow the measurements of the 
hyperfine splitting of electronic states in atoms with an 
accuracy of the order 10~ 6 , that provides the possibility 
of hyperfine anomaly studies. 

The hyperfine splitting is sensitive to the magnetic cur- 
rent in the nucleus, and the hyperfine structure (hfs) con- 
stants extracted from the experimental data are related 
to the matter distribution of the nucleus. Therefore, the 
measurements of the hfs anomaly in neutron-rich halo 
nuclei can give a unique way to investigate the neutron 
distribution and the cluster structure. 

For the experimental values of the magnetic moment 
of the Be isotopes, we refer to the measurements in Refs. 
0. Q Q But still there are no experimental data on the 
hyperfine anomaly for these nuclei. 

A theoretical study has been performed by Fujita et 
al. in order to verify how the halo structure man- 
ifests itself in the hyperfine structure. They calculated 
the hfs anomaly for Be isotopes both in the core+neutron 
model and in the single particle model, where the mag- 
netic moment is well reproduced. They found that the hfs 
anomaly for 11 Be is large compared to that for the 7 ' 9 Be 
isotopes, and that this was indicative of an extended neu- 
tron distribution and a halo structure in 11 Be. 
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However, the approach used by Fujita et al. in Ref. 
suffered from a poor knowledge of the ground state 
wave function of the 9 ' n Be isotopes. The value of the 
hfs anomaly e defined by the folding of the electronic and 
nuclear wave functions is rather sensitive to the spatial 
distribution of the valence neutron wave function and to 
the weights of the possible mixed states in the description 
of the ground state wave function. During the last decade 
new information on the 11 Be ground state wave function 
became available from the cross section measurements in 
the p( 11 Be, 10 Be)d reaction and the weight of the 2si/2 
state admixture was found at the value 0.84 Q which is 
close to the theoretical estimation. 

Besides this, the cluster (core+neutron) model while 
being rather good for the 11 Be nucleus, might fail for 
9 Be, which has essentially a three-cluster (a + a + n) 
structure On the other hand, there is no strong ev- 
idence for the three-body structure of 9 Be in its ground 
state. So this is another interesting problem which can 
also be investigated in hfs anomaly studies. 

The n Be 1/2+ ground state can be represented by an 
admixture of the 2sx/2 and ld 5 / 2 valence neutron states, 
w Be(0+)®v{2s 1/2 ) and 10 Be(2+)®i/(ld 5/2 ). The weight 
of the s- state, related to the spectroscopic factor of this 
configurations obtained in both the shell model and the 
simple excited core cluster model, are found between 0.5 
and 0.8 (for more information on the spectroscopic fac- 
tors we refer to Refs. j^ElE!)- In our calculations we 
take the weights of the states w(2s 1 / 2 <8 + ) = 0.72 and 
w(ldt.i') 2 + ) = 0.28, consistent with the experimental 
data |12| (for more details, see 0]). 

Let us mention that in Ref. 6] the hfs anomaly 
has been calculated within the Bohr-Weisskopf approach, 
where the electronic wave function is found for the in- 
terior of the nucleus. It is not applicable to the 11 Be 
nucleus, where a significant part (88%) of the valence 
nucleon wave function is outside of the range of nuclear 
potential, indicating a halo ^ n particular, the 10 Be 
core root mean square (rms) radius 2.61 fm is rather 
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small compared to the mean core-neutron distance (6 
fan). Thus, in Ref. Q a rather simplified electronic wave 
function is used which significantly differs from the di- 
rect numerical solution of the Dirac equation in the re- 
gion r < 10 fan. Besides this, in calculations of the elec- 
tronic wave function, the electron screening effects of the 
Coulomb potential were not taken into account, although 
this effect can also be essential in the calculations of the 
hfs anomaly. 

In the present paper, the relativistic electronic wave 
functions are calculated for the extended nuclear charge 
density distribution defined by the cluster structure of 
the nucleus [lj|. We assume that two electrons in the Is 
state form the closed shell and are relatively unperturbed 
by the third electron in the 2s state, and we define the 
wave function of this third electron as a solution of the 
Dirac equation for a nuclear charge potential screened by 
the Is 2 electronic closed shell [lfij . 

We compare our electronic wave function with that 
of Fujita et al. 6]. We analyze the hfs of 9 > n Be with 
regard to the available experimental information using 
more realistic descriptions of both the nuclear and the 
electronic wave functions. The aim of this paper is to 
answer the questions whether the hfs anomaly in these 
isotopes reflects the halo structure, and whether the hfs 
anomaly for n Be isotope is larger than for 9 Be. 



and a e is the Dirac matrix for relativistic electrons. 

For a .ZV-electron system, the hyperfine interaction 
Hamiltonian can be defined as 



H = — 



JW-j(r') 



d 3 r d 3 r'. 



(4) 



With the Neumann expansion of | r i; r / | m equation J3J , 
the interaction Hamiltonian takes the form: 



H = -~ / /VH(O[T(r<)-U(r>)] (rr0 dVrf 

\ J JWoeFr^j-uy^A (5) 



where 



T(r<) • U(r>) = E E T ^( r <) t/ A,(r>) (6) 

A v 



A v 



and 



II. MAGNETIC HYPERFINE STRUCTURE 

The magnetic hyperfine interaction Hamiltonian is de- 
fined by 



H = - J J(r)-A(r) d 3 r. 



(1) 



Here, J is the nuclear current density, A is the vector 
potential created by the atomic electrons 



A(r) = \ J 



1 r ^ ^ 



(2) 



where j(r') denotes the electron current density operator 



j(r') = - e a e 5{r - r') 



(3) 



4tt 



2A + 1 



Y Xv {f) 



with r < and r> being the smallest and largest value of the 
nuclear (r) or electronic (r e ) coordinates. The A = 1 term 
is the magnetic dipole interaction between the magnetic 
field generated by the electrons and the nuclear magnetic 
dipole moment due to the extended nuclear matter dis- 
tribution. The A = 2 term is the electric quadrupole 
interaction between the electric field gradient from the 
electrons and the non-spherical charge distribution of the 
nucleus. 

The hyperfine interaction couples the electronic angu- 
lar momentum J and the nuclear one I to an hyperfine 
momentum F = J + I. The magnetic hyperfine splitting 
energy W for a state | IJFMp = F > is defined as the 
matrix element of the Hamiltonian Ti, 



W {IJ)FF = < IJFF | U | IJFF > 

= < IJFF | IM, Jm>< IM, Jm\H\ IM', Jm! >< IM', Jm! \ IJFF > 



(7) 



m.nil 



with the matrix element 



< IM, Jm\H\ IM', Jm' >= -2i- ^ < IM | / J(r) ■ (A° x (r) + A c x {r))d 3 r | IM' >, 



(8) 
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where 



V 

is the expression for a point nucleus, and 



dr e 



(re) 



A+l 



F Kj (r e ) G Kj (r e ) < Jm \ C\(f e )u e \ Jm' > 



Ai(r) = 



(-yC*(f) I dr e F Kj (r e )G Kj (r e ) ( 



IKj { 



„A+2 
_e 

„A+2 



„A-1 



—)< Jm\ Ct v {r e ) er e \ Jm 1 > 



is the correction term due to the finite extension of the 
nuclear density. 

The functions F kJ , G kJ are the radial parts of the 
large and small components of the Dirac wave function 
of the electron, with the quantum number n = ±( J + i) 
for J = L =p i and the orbital angular momentum L. 

In the dipole approach (A = 1) the expression for 
W(u)ff reduces to 



W [IJ)FF = < IJFF\H\IJFF > 



(9) 



l[F(F + l)-I(I + l)-J(J + l)] ai 



where aj is defined by 
ai = t^t- — r < II I ? 



A 



{M l z {v l ) + M s z {v l ))\II> 



with the Z components of the angular (or spin) magnetic 
moments M^ s )(rj). Here, the summation runs over all 
nucleons. 

Taking into account the density expansion of the nu- 
cleus, these moments are 



M ( (r 4 ) = g\U[ 



M s (r,) = gi[ Si 



F kJ G kJ dr 



F Kj G Kj dr 



F kJ G kJ (-) 3 dr], 

n 



D, 



F Kj G Kj (-fdr], 

Ti 



where Dj = - J § [s 1 <g> C 2 (fi)} . 

Thus ai can be expressed through the hfs constant for 
a point nucleus as 



(10) 



where e is defined as the hfs anomaly in the Bohr- 
Weisskopf effect and 6 is the Breit-Rosenthal-Crawford- 
Shawlow ('BRCS') correction 

The hfs constant for the point nucleus is 



(o) 2en fiN/i 



/J(J + 1) 



F « J (r)G% J (r)dr, (11) 



where (i =< II | ^ g l s Si + g\li \ II > defines the mag- 

i=l 

netic moment of the point nucleus in nuclear magneton 
units hn- The functions i 7, Kj (r) , Gq (r) are the ra- 
dial parts of the large and small components of the Dirac 
wave function for the electron in the point nucleus ap- 
proximation. 

So, the hfs anomaly is defined by 



e = -~ 1 £{<II\(9? ) Si + 9?k)K a (ri)\II> 
M i=i 

< II\(gi l) U + D t )K b (r t ) | II >} . (12) 

Here, b = [/ °° F£ j Gq J dr]~ 1 is a constant obtained from 
EH and 



K a (n) 



• i 

1 



F Kj G Kj dr 



F Kj G Kj (-fdr. 
r, 



(13) 
(14) 



The 'BRCS' correction is 8 = 1 - bK a (oo) [T^. 

Let us assume the two-cluster nuclear wave function 
defined as a superposition of different configurations 
$[ s = [^i" s ® ^J 7 , associated to the coupling of a 

core state <J> ; c and the valence particle wave function 

ifp, . J c , s c and j, s are the total angular momentum 
and spin of the core fragment and valence nucleon, I is the 
orbital angular momentum of their relative motion. For 
each configuration, the contribution to the hfs anomaly 
is given by 

« - ~Y,[< II \^s i+9 ^l)K^R)\II> 



M 



M 



- <II\(gl i) ^l + D t )K b ( 1 ^R)\II> 



M 



M 



(15) 



where indices correspond to the core (i — 1) and valence 
nucleon (i = 2), M = m\ + m-i is the mass of the whole 
system; r and R define, respectively, the electronic and 
relative radial coordinates of the two nuclear fragments. 
gjp and g^ are the gyromagnetic ratios of the i-th frag- 
ment orbital motion and spin, respectively. 
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III. ELECTRONIC WAVE FUNCTIONS 

In the experiment pj, the hyperfine splitting of the 
electronic levels (F) is measured for the ground state of 
Be + ions. The electronic ground state of the Be + ion is 
represented by the ls 2 2s configuration, and the hyperfine 
splitting of the states F — 2 and F = 1 is found at about 
1.256 GHz. 

In the present paper, we calculate the hfs anomaly 
of the Be isotopes using the electronic wave functions 
pnj _ /W ^ q K j _ sW obtained as solutions of the 

r r 

Dirac equation [18j taking into account the (Is 2 ) electron 
screening and the extended density of the nuclear charge 
distribution. The electronic wave functions in the region 
r < 22 fm are obtained for the Coulomb potential defined 
by the expression 



where p e i (x) , is the electron density distribution of the 
two Is electrons in the closed shell. This potential can 
be approximated by 



V(r) 



Zahc 1 



te- 



ll 



(19) 



where k is fitted to reproduce the initial potential. So, 
the effective charge for the 2s electron in the Be + ion is 



J eff 



1) 



(20) 



with x = 0.000068 fm' 1 . 

Finally, the electronic wave functions calculated in the 
interior of the nucleus are matched to those obtained with 
the screened Coulomb potential. 



V(r) = 4tt 



Zahc 



px 2 dx + r / p x dx > , (16) 



where p(x) is the charge distribution of the nucleus. This 
distribution is obtained in the cluster model of the Be iso- 
topes using the ground state nuclear wave functions. To 
calculate the electronic wave functions we use the numer- 
ical methods suggested in Ref. [l9| . 

The asymptotic electronic wave functions (r > 22 fm) 
arc found in the form 



f(r) 
9(r) 



-Dr 



N 

E< 



-Dr 



aZ, 



eff 



E^ s+ " 



where A = ±(E N + mc 2 ), B = ^(-E N + mc 2 ), D = 
\J~AB, s = yjk 2 — (aZ e ff) 2 and the 2s electronic energy 
En given by 



Ej\[ = mc 



Z eff a 



N + ^/k 2 -(Z eff a) 2 



-1/2 



(17) 



with N = 0,1, 2..., and k = ±1, ±2, ±3..., are found with 
the effective charge Z e ff depending on radius r. a v and 
b v are expansion coefficients (see Ref. [T^l- 

In the case of a Be + ion, there are three electrons or- 
biting the nucleus. We can make a simplification by con- 
sidering the 2s electron in the Coulomb potential of the 
nucleus screened by the closed electronic shell Is 2 . Thus, 
the asymptotic electron wave functions are defined for 
the screened Coulomb potential written in the form: 



V(r) 



Zahc 



1 — 2tt < / p e ix dx + r I p e ixdx 



+ l -{U/^pT{r) 



(18) 



IV. RESULTS AND DISCUSSION 

The hfs anomaly e for the Be + ion is obtained from 
the relation aj = 0/(1 + e + 5) where is the Fermi- 
contact parameter found for a point nucleus |20j . In the 
corc+nucleon model of the Be isotopes the values aj and 
e are determined by the valence nucleon wave function 
(that, in general, is an admixture of different single parti- 
cle states) and the electronic wave function depending on 
the charge distribution of the nucleus. Hence, the amount 
of clustering or the existence of a "halo" in the Be iso- 
topes nuclei can be analyzed through the "hfs" anomaly 
evaluation. 



We consider the 11 Be nucleus as a two-body system 
composed of a 10 Be core nucleus in different states, + 
and 2 + , and a valence neutron. The 11 Be ground state is 
described by a superposition of the neutron states, 2s!/ 2 
and ld 5 /2i as 



= /3|[ 10 i?e(0+)®7^ 1/2 ] 1/2+ ) (21) 
+ \\[ w Be(2 + )®n d5/2 ] 1/2+ ) 

with the weights w s = /3 2 and Wd = A 2 of the s- and d- 
waves obtained in [!2 , which allow a good description of 
the 11 Be interaction and break-up reaction cross sections 
at intermediate and high energies (see Ref. ^3)- 

In our approach, the magnetic moment of n Be in the 
ground state is 



/' 



W s p s 

[ I72 ? 

W a p s 



E E 

m B .mi nij ,Ms 



r<2 2 



M (2+) 



M s 
7 



sm s lmi SMsjmj 



mi. 



(22) 
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WdUd + TrWdPc - -W d p(2 + ), 

15 6 
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TABLE I: Core and valence neutron parameters used for the 
calculation of the core-nucleon wave function and core density 
parameters. E x is the excitation energy of the core state, Wj 
are the weights of the different configurations , Vo, Ro are 
the parameters of the Woods-Saxon potential with diffuseness 
parameter ao = 0.5 fm, B n is the neutron separation energy. 



E x 


core 


neutron 




Woods-Saxon potential 6 




Jc 


nlj 


Wj 


Vo Ro 


B n 


(MeV) 








(MeV) (fm) 


(MeV) 





0+ 


2Si/ 2 


0.7f9 


-70.83 2.483 


0.504 


3.368 


2+ 


1^5/2 


0.281 


-89.22 2.483 


3.872 



6 The 10 Be charge rms radius is related to the matter radius 
r m as rms^x/r^ + 0.64, where r m is taken from Ref. |23| 



where (s, m s ), (I, mi), and (j, rrij) are the valence neutron 
spin, orbital angular momentum and total momentum, 
and their projections. S, Ms are the spin of the 10 Be 
core nucleus and its projection. The magnetic moment 
of the 10 Be core fragment related to its orbital motion is 
P c = A(I-i) ^ A = 11. Here, /i s , p, d , and fi(2 + ) denote 
the magnetic moment of the neutron in the Si/ 2 and d 5 / 2 
states and the 10 Be core in the 2 + state. We use the 
values fi s = [id = 5.9s"' = — 1.9135/ijv- The magnetic 
moment of 10 Be in the excited state V = 2 + has been 
calculated in Ref. |2jJ in the shell model, and has been 
found to be A*(2 + ) = 1.787/ijv- 

The ground state wave function (|21|l gives the 11 Be 
magnetic moment value — 1.784/ijv that is larger than 
the experimental value /ij = — 1.6816(8)/zat |5j. This 
problem with the magnetic moment description was dis- 
cussed, in particular, in and the weights of the s- 
and d- waves have been found w s =0.55, and 11^=0.45. 
In our simplified model, we reproduce the experimental 
value of the 11 Be magnetic moment with u> s =0.5, and 
uu=0.5, which are the same as those used in calculations 
|6j . This example gives a measure of the sensitivity and 
the precision of the calculations. 

The valence neutron wave function for each partial 
state in 11 Be is obtained as a solution of the Schrodinger 
equation with the same Woods-Saxon potential parame- 
ters as in Ref. , and reproduce the neutron separation 
energy B n . These parameters are given in Table [|] 

The 10 Be core density is parametrized using the har- 
monic oscillator model [15J 



P( r ) = Po(l + l{r/a) 2 )exp{-(r/a) 2 ), 



(23) 



where po is a normalization factor, a — 1.856 fm, and 
7 = 0.610. 

The charge density distribution of the 11 Be nucleus is 
defined as 



p(r) = Po 



J d 3 x Pc(|r-_Lx|), 



where po is a normalization factor, p c is the core density 
and §l( x ) defines the radial part of the wave function for 
the valence neutron. The charge density distribution in 



(|16fl and the electronic density distribution in (|18fl give 
the Coulomb potential entering the Dirac equation for 
the wave function of the outer 2s electron. 

With the weights from Table [I] and parameters men- 
tioned above we obtain the value of the hfs anomaly for 
the 11 Be + ions as 



b 

Pi 



>M.{AS(§)-(1-Ei«^)} 

~M 

: a(™)+£(2+)jft(™)} 

iy M' y ' lK M Jf 



^A»M^)-(l-=i)^(§)} 



(24) 



where 



- 1a-m*5(§)-aS(§)}. 



1/2 



3 3 
±-- 

4 j + 1 



for j = I ± |. The value of S(2+) is taken as S(2+) = 
— 1.034 (see Ref. Here, m and M c are the mass of the 
neutron and of the core 10 Be, respectively. The ratio ^jf- 
(jt) takes into account the motion of the core (valence 
neutron) relative to the 11 Be center of mass (cm). 
The terms Kf and K\ are expressed as 

oo 

K?( mi ) = J \<5>i{R)\ 2 R 2 dRK a {m t R) (25) 
o 

oo 

K?{mi) = J \>5>i{R)\ 2 R 2 dRK h {m l R) (26) 
o 

where K a and K b are defined in |T3J and ljl"4"|) and &i(R) 
is the radial part of the valence nucleon wave function, 
for the state with the orbital angular momentum I. 

Notice, that the terms and K\ in l|24|l are calcu- 
lated for the p - wave nucleons composing the excited 
10 Be nucleus and contributing to the magnetic moment 
of the excited 2 + state. In these calculations we used the 
folding of the electronic part K a ^ with the 10 Be den- 
sity distribution, assuming that the spatial distribution 
of the p - wave nucleons is similar to the charge density 
distribution in 10 Be. 

Table |H] shows the values of the hfs anomaly calcu- 
lated without screening effect (e" s ), or that for a homo- 
geneously charged sphere with the rms radius 2.61 fm 
(e hs ), and the value e obtained by including all the ef- 
fects. The contribution of the terms Kf and K\ in l|24() 
to the hfs anomaly are given in the row named " core" . In 
the Table, r n is the mean distance of the valence neutron 
from the 11 Be cm. 

One can see from the Table, that the hfs anomaly is 
sensitive to the screening effect (about 25%). The total 
value weakly depends on the shape of the charge density 
distribution (about 3 %), but the effect is different for 
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TABLE II: The hfs anomaly in calculated for different 

states of the 10 Be core. B n is the neutron separation en- 
ergy, r n is the neutron distance, e na is the hfs anomaly value 
obtained without consideration of the screening, e is calcu- 
lated for a homogeneous spherical charge distribution, and e 
is obtained including all the effects. In the last two rows the 
weighted sums of the values are given. 



0.03 



State 




ns 


hs 




r n c 


e c 


nlj 


(MeV) (fin) 


(%) 


(%) 


(%) 


(fin) 


(%) 


2S1/2 


0.504 6.056 


-0.0558 


-0.0694 


-0.0719 


6.165 


-0.120 




3.872 2.634 


-0.0090 


-0.0103 


-0.0112 


3.551 


-0.0233 


core 


2.623 


-0.0006 


-0.0007 


-0.0008 


2.588 




Total, 


w s = 0.72 


-0.0433 


-0.0534 


-0.0556 






Total, 


uis = 0.5 


-0.0335 


-0.0410 


-0.0430 




-0.0717 



cited from Ref. 



distinct partial states, and is more pronounced for the d 
- wave neutron (about 8%). 

These values of e can be compared with those from 
@, e(2s 1/2 ) = -0.120% and e(ld 5/2 ) = -0.0233%. The 
values from 6] are obtained without screening effect in a 
simplified model of the nuclear charge distribution, ap- 
proximated by a homogeneously charged sphere with the 
radius 1.2A* = 2.669 fm. 

The value of the 'BRCS' correction is S = 0.0476% 
that is comparable with e. 

The hfs anomaly value for each partial wave function 
strongly correlates with the neutron separation energies 
B n and r n . Thus, the hfs anomaly value, obtained for the 
s-wave neutron is larger than that for the cZ-wave by an 
order of magnitude. Hence, the value of the hfs anomaly 
is very sensitive to the weights of the partial waves. For 
example, with the weights from Ref. w s = Wd = 
0.5 we get the value e = -0.0430% which is essentially 
smaller than that obtained in Ref. @, e = -0.0717%. 
The contribution of the core fragment is also small. 

Thus, the precise measurement of the hfs anomaly 
might shed light upon the spectroscopic factors of the 
s- and d- states and the origin of the i + state of 11 Be. 

The halo nucleus charge distribution might signifi- 
cantly deviate from a spherical distribution, defining the 
radial dependence of the electronic wave function. At 
the same time, the 10 Be— n mean distance is 6.7 fm, so 
that r n ~ 6 fm, a value exceeding a few times the 10 Be 
core rms radius. Thus a significant part (88%) of the 
valence nucleon wave function is outside of the range of 
the nuclear potential. Hence, the result of the folding 
of the electronic and nucleon wave functions essentially 
depends on the radial behavior of the wave functions. 

Fig n shows the radial dependence of the electronic 
wave functions inside the 11 Be nucleus, calculated in the 
Bohr-Weisskopf approach 0, the simplified approach 
, and that obtained here as a numerical solution of the 
Dirac equation. One can see that the calculated wave 
functions are close to each other at r < 3 fm and sig- 
nificantly different at higher radii. For 11 Be, the folding 
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0.02- 



0.01 



0.00 



-num 
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F[6] 




r(fm) 



FIG. 1: Electronic wave functions (J=|,/t = — 1) inside the 
11 Be nucleus, obtained in the Bohr-Weisskopf approach |24| 
(BW - dashed line), the simplified approach of Fujita et cd. 
[g (F- dotted line), and the numerical solution of the Dirac 
equation (num - solid line). 



of the electronic wave function with the wave function of 
the weakly bound valence nucleon gives significantly dif- 
ferent results. Clearly, with the wave function of Fujita 
et al. [6|, the hfs anomaly is overestimated for weakly 
bound nuclei. 

This overestimation, in particular, can be seen from 
comparison of the hfs anomaly values for the s-wave va- 
lence neutron with close r n values. The result of our 
calculations is smaller than that from Ref. [jj]. For the 
d-wave neutron, with different r„ values we obtain the 
hfs anomaly value twice as small as that from Ref. §. 

Besides the difference in the electronic wave function 
radial dependence, the assessment [24[ of the integrals 
with the electronic wave function in (|25|l and 126|) as 
- KUtt) = CK?{r$) with C = 0.68 is rather 
crude, and the coefficient C for different partial waves 
varies within 20%. 

To compare the calculated value of the hfs anomaly in 
11 Be with that for 9 Be, as the next step of our analysis, 
we perform similar calculations for the 9 Be + ion. 



J Be 



The hfs anomaly for 9 Be can be calculated within the 
two-cluster model of 9 Be. This nucleus can be regarded 
as a system composed of a p3/2-wave neutron and a 8 Be 
core in the ground (0 + ) and excited (2 + ) states. 

Note, that in general, the states + , 2 + , and 4 + are 
realized in 8 Be, but only + and 2 + contribute to the hfs 
anomaly in 9 Be (corresponding to the valence neutron 
angular momentum 1 = 1). The contributions with I = 3 
are expected to be small. 

The 9 Be ground state is given by a superposition of 
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the + and 2 + states 



= p\[ 8 Be(0 + )®n P3/2 } 3/2 -) (27) 
+ A|[ 8 Be(2+)0n P3/2 ] 3/2 -). 



With this wave function, denoting u> + = (3 
w 2 + — A 2 , the magnetic moment of 9 Be 
core+cluster model is 



and 
in the 



r 1 3 

// = w +{fi p + fi c } + w 2 +{-fi p + -fj, c 



1 



M(2+)}. (28) 



Here, indices in ujjjt correspond to the core states — 
0+,2+. 

The magnetic moment of the 8 Be core related to its 
orbital motion is [i c = ma-i) ^" 1 ^ = 9. Here, fi p , and 
/x(2 + ) denote the magnetic moment of the p 3 / 2 neutron 
and the 8 Be core in the 2 + state, pL v — — 1.9135/ijv an d 
M(2+) = 1. 

With the weights w + — 0.535 and w 2 + = 0.465 ob- 
tained with the spectroscopic factors from Ref. |2^| the 
magnetic moment is fi = —1.0687 /x;v- The change of the 
parameters to w + — 0.579 and w 2 + — 0.421 allows one 
to reproduce the experimental value of the 9 Be magnetic 
moment fj,i = — 1.1447/ijv. This again gives a measure of 
the sensitivity of the calculated value to the weights of 
the states. 

In our approach, these two states of 9 Be with the core 
in the ground and excited state are characterized by dif- 
ferent neutron separation energy B n listed in Table ITTll 
For the case with 9 Be (0 + ), the valence neutron wave 
function is obtained as a solution of the Schrodinger equa- 
tion for the Woods-Saxon potential Vb = —43.61 MeV, 
a = 0.5 fm, and R = 2.46 fm giving B n = 1.665 MeV 
and the 9 Be charge rms radius 2.519 fm. 

In the case of excited 8 Be in the resonance state 2 + , 
the parameters of the Woods-Saxon potential are taken 
as Vb = -49.80 MeV, a = 0.5 fm, and R = 2.49 fm to 
reproduce B n (see Table lITTll . 

The 8 Be core density is parametrized in the harmonic 
oscillator model l|23(l . with the parameters a — 1.749 fm 
and 7 = 0.619 for the + state, and a — 1.768 fm and 
7 = 0.624 for the 2 + state, respectively. The density 
parameters give the rms radius of the core equal to Rq 
for each state. 

The hfs anomaly for 9 Be + ion is 



+ 



(w 0+ + ^w 2 +)[i p {K^(^) 

|A W2+{M (2+)i? f (^)+E(2^ 
5 [ii M 



-(1-E*)*f(^» 



TABLE III: The hfs anomaly e in 9 Be+ calculated for dif- 
ferent states of the 8 Be core here and in (fj, wj£ are the 
weights of the different configurations. B n is the neutron sep- 
aration energy, r n is the neutron distance. In the last row the 
weighted sum of the values is presented. 



)KH — )} 



The values Kf and K\ are calculated with (|23|l and (|2*r]|l . 

The values Kj" and K\ are calculated with the density 



Core 






r n 




r„ a 


e a 


Jo 


(MeV) 


(fm) 


(%) 


(fm) 


(%) 


0+ 


0.535 


1.665 


3.200 


-0.0440 


2.569 


-0.0249 


2+ 


0.465 


4.705 


2.630 


-0.0066 






core 


0.535 




2.490 


0.0063 






Total 








-0.0236 




-0.0249 



a cited from Ref. 



distribution of the core fragment in the 9 Be nucleus, tak- 
ing into account the 8 Be-neutron relative motion. 

The calculations show that the contribution of the 
term T,(2 + )K is small and we can neglect it. 

In Table ITTT1 the hfs anomaly values are listed for each 
8 Be state. The total value, obtained as the weighted sum 
is given in the last row of the Table. It is close to the 
values e = -0.0249% @, e = -0.0243% 0. 

The value of the 'BRCS' correction is S = 0.0451% 
that is larger than e. 

One can see that the result for 9 Be is twice as small 
as that for 11 Be. This corresponds to the conclusion in 
Ref. 0, that the value of the hfs anomaly reflects the 
extended neutron distribution in 11 Be and might indicate 
a neutron halo. At the same time, the difference in the 
hfs anomaly values is not as large as in Ref. [(|. 

As with 11 Be case, the 9 Be hfs anomaly values for each 
partial wave function strongly correlate with the neu- 
tron distance r n from the Be cm. In particular, the hfs 
anomaly obtained for the p-wave neutron for the + state 
of the 8 Be core is larger than that for the 2 + state by an 
order of magnitude. The contribution of the core frag- 
ment is also small. 

We should note, that we use a rather simplified model 
of 9 Be which most likely has a three-cluster structure 
(a + a + n). In the three-cluster approach, there are 
partial states, contributing to the 9 Be ground state wave 
function, which are not taken into account in the two- 
cluster model, and which also contribute to the hfs 
anomaly value. 



CONCLUSION 



In the present paper, we have calculated the hfs 
anomaly in the 9,11 Be isotopes. The nuclear wave func- 
(29) tions have been calculated in the simplified corc+ncutron 
model. In calculations of the realistic electronic wave 
function, the charge distributions of the clustered nu- 
cleus and electrons in the ls 2 2s configuration in Be + are 
taken into account. 

It is found that the value of the hfs anomaly for 11 Be 
(e = -0.0556% for w s =0.72 and e = -0.0410% for 
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w s =0.5) is larger than that for 9 Be (e = -0.0235%), 
that agree with the conclusion of Fujita et al. |{J. 

For 9 Be we reproduce the hfs anomaly value, obtained 
in [(|, e = —0.0249%, and for 11 Be we get the values with 
different weights Ws of the s-wave, which are essentially 
less than that in |g|, e = —0.0717%, obtained with the 
weight w s =0.5. 

In calculations of the hfs anomaly we did not con- 
sider the three-cluster structure of 9 Be. We did not take 
into account the probable contribution of other states ad- 
mixed in the ground state wave function of 9 Be, which 



also might influence the calculated value of the anomaly. 
Thus, the difference in the hfs anomaly in the 9 Be and 
11 Be isotopes might be even less. 

Finally, we can conclude, that, on one hand, the hfs 
anomaly correlates with the neutron spatial distribution. 
As the anomaly value is sensitive to the ground state 
wave function of the nucleus, there is no unambiguous 
correlation between the hfs anomaly value and valence 
neutron distribution, and each nucleus requires a sepa- 
rate investigation. 
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